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CNJ ' Abstract. A GRMHD model of disk outflows with neutrino-driven mass ejection is presented. The 

model is used to calculate the structure of the outflow in the sub-slow magnetosonic region and 
the mass loading of the outflow, under conditions anticipated in the central engines of gamma-ray 
bursts. It is concluded that magnetic launching of ultra-relativistic polar outflows is in principle 
possible along low inclination field lines (with respect to the symmetry axis), provided the neutrino 
luminosity is sufficiently low, L v £ 10 52 erg s _1 . 
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INTRODUCTION 



> 

Neutrino-driven mass loss plays an important role in the late phases of proto-neutron 
star evolution 11E-B3], and in the outflows expelled from the hot, dense disks in hyper- 



accreting black hole systems (e.g., Refs. |al2D)- In weakly magnetized systems those 



neutrino-driven outflows are subrelativistic, because the entropy per baryon generated 



by neutrino absorption at the base of the wind (s/kg ^ 100) does not reach the values 



required to accelerate the outflow to high Lorentz factors. However, in highly magne- 



tized systems the rotational energy of the system can be extracted magnetically and the 
outflow can in principle accelerate to high Lorentz factors provided the mass loading of 
magnetic field lines is sufficiently low [ 

Magnetic extraction of the rotational energy of accreted matter in collapsars or, al- 
ternatively, a highly magnetized proto-neutron star is a plausible production mechanism 
for the y-ray emitting jets inferred in long GRBs. To account for the characteristic lumi- 



nosities and Lorentz factors observed in those systems two conditions must be fulfilled: 
(i) an ordered magnetic field with strength of the order of 10 14 — 10 15 G must be present 
at the flow injection point and (ii) the ratio of the neutrino-driven mass flux to magnetic 
flux must be sufficiently low. The possible existence of such strong magnetic fields in 
systems like those mentioned above appears to be indicated by observations, and is also 
supported by recent numerical simulations. However, the dependence of the neutrino- 
driven mass loading on the conditions in the disk and on the geometry of the magnetic 
field is yet an open issue. 

A common approach in studies of magnetized disk outflows is to seek self-similar 
solutions of the trans-field equation (e.g., Refs. USEilEI)- While this approach allows 
the magnetic field geometry to be calculated self-consistently, incorporation of gravity 
(in the relativistic case) as well as neutrino heating in the flow equations is precluded. 
Furthermore, it does not allow simple matching of the self- similar outflow solution 
to a Keplerian disk (but c.f., Ref. Ill 211 \ The immediate implication is that the mass- 



to-magnetic flux ratio, which is determined by the regularity condition at the slow 
magnetosonic point, remains a free parameter. One must therefore adopt a different 
approach. One method is to seek solutions of the flow equations that pass smoothly 
through the slow magnetosonic point, assuming the magnetic field geometry is given. 
Such an analysis has been carried out recently for winds from a disk around a hyper- 
accreting black hole [9]. A brief account of the model and results is given below. 



THE MODEL 

We consider a stationary, axisymmetric MHD wind expelled from the surface of a 
hot, magnetized disk surrounding a non-rotating black hole. The range of conditions 
in the disk is envisioned to be similar to that computed by Popham et. al. 111411 for 
hyperaccreting black holes with mass accretion rates 10 -1 — 10 M s _1 and viscosity 
parameters oc v ,- 5 = 0.01 — 0. 1 . Under such conditions the dominant cooling mechanism in 
the disk is neutrino emission. The major fraction of the neutrino luminosity is generated 
in the inner disk regions, within 10 Schwarzschild radii or so, and lies in the range 
L v = 10 51 - 10 54 erg s _1 for the above range of accretion rates and viscosity parameter. 

The model calculates the structure of the GRMHD outflow below the slow magne- 
tosonic point for a given magnetic field geometry, treating the neutrinos emitted from the 
disk as an external energy and momentum source. To simplify the analysis the neutrino 
source is taken to be spherical with radius R v . The model is characterized by three pa- 
rameters: the black hole mass Mbh, the neutrino luminosity L v , and the neutrino spheric 
radius R v . It solves a set of coupled ODEs that are derived from the general relativistic 
energy-momentum equations, describing the change along a given streamline *P(r, 6) — 
const of the specific energy $ , specific entropy s, and poloidal flow velocity u p : 

(p/m N )k B Ts' = -u a q a , (1) 
pcV = -q u (2) 

(\nu p )' = ^. (3) 

Here ' denotes derivative along streamlines u a d a , q$ denotes the source terms associ- 
ated with energy and momentum exchange with the external neutrino source, u a is the 
outflow 4-velocity, p is the baryon rest mas density and T is the temperature. The de- 
nominator on the R.H.S. of Eq. © is given explicitly as D = —(a 2 — R 2 £l 2 —M 2 ) 2 (u 2 — 

u sm) ( m p — u fm) I M l' where ua, usm and ufm are the Alfven, slow and fast magnetosonic 
wave speeds, respectively, a is the lapse function, Q. is the angular velocity defined be- 
low, R is the cylindrical radius and M = u p /ua is the Alfven Mach number. The term F 
can be expressed as F = ^QnBp)' + £ 2 (lna)' + £ 3 (ln/2)' + £ 4 (ln#)'+ £ 5 (]ns)', where 
the coefficients ^ are functions of the flow parameters, given explicitly in Ref. 0, 
and B p is the poloidal field component. Since the derivatives (lnB^)' and (\nR)' de- 
pend on the magnetic field geometry which is unknown a priori, additional equation is 
needed. Our approach is to invoke a given field geometry. To examine the dependence 
of mass flux on the latter, we obtained solutions for different magnetic field configura- 
tions, focusing particularly on split monopole and r self-similar geometries. Equations 



©-© are augmented by an equation of state for the mixed fluid of baryons, photons and 
electron-positron pairs. In addition there are three integrals of motion of the MHD sys- 
tem: the mass-to-magnetic flux ratio T7 (*F), the angular velocity of magnetic field lines 

and the specific angular momentum Jzf The two invariants and Jzf 
are fixed by a choice of boundary conditions. The mass flux 77 (^P) is an eigenvalue of 
the problem. The integration starts sufficiently deep in the disk atmosphere where the 
density is high, the flow velocity is very small, and the ratio of baryonic pressure and 
the light fluid pressure is about unity so that the specific entropy is dominated by the 
baryons. The value of r\ is adjusted iteratively by repeating the integration, changing the 
boundary value of the poloidal velocity u p q each time, until a smooth transition across 
the slow magnetosonic point is obtained. 



RESULTS 

Typical solutions are shown in Figs 1 and 2. Those were obtained using a split monopole 
geometry (see Ref. |9] for more details). In Fig. 1 the change along a field line of the flow 
quantities indicated is plotted against the normalized height above the disk midplane, 
z/Rq, where Rq is the radius at which the field line meets the disk. The inclination angle 
of the field line in this example is about 12° with respect to the rotation axis, so it lies in 
the so-called stable regime. The mass flux is thermally driven in this regime rather than 
centrifugally driven. As seen, the slow magnetosonic point is located at z — Rq, which we 
find to be typical for injection along low inclination field lines. The integration started 
well beneath the slow magnetosonic point where the specific entropy is dominated by 
the baryons. In this region adiabatic cooling is negligible and kinetic equilibrium is 
established. The net energy deposition rate nearly vanishes, viz., q'(To) ~ 0, and the 
temperature is Tq ~ T v , where T v is the neutrino equilibrium temperature. As the flow 
accelerates its temperature starts falling, and since the neutrino cooling rate is very 
sensitive to the temperature q l increases rapidly leading to the steep rise of the entropy 
per baryon seen in Fig 1. The terminal value of s is however modest, s/kg < 100 in all 
cases examined, so the net outflow energy is practically not affected. The main effect 
of neutrino heating is not to change the net outflow energy along streamlines, but rather 
to enhance the baryon load by increasing the light fluid pressure in layers of higher 
baryon density. Thus, the initial value of <S presents, essentially, an upper limit for the 
asymptotic Lorentz factor of the wind. The values of the mass flux, which we somewhat 
arbitrarily define as M = Pqu p q1tzR^ (for a two-sided outflow), that corresponds to the 
cases depicted in Fig. 1 are 10 28 and 7 x 10 30 g s _1 for T v = 2 and 4 MeV, respectively, 
and the the corresponding values of the total energy per baryon are <§ = 5 x \0 2 m p c 2 
and 2m p c 2 for a surface poloidal field strength of 10 15 Gauss. 

The dependence of the mass flux on the inclination angle of the field line is shown 
in Fig. 2. As seen, the mass flux is quite sensitive to the inclination of the field line. 
For field lines in the unstable regime (inclination angles larger than 30° to the vertical) 
the wind is centrifugally driven and, for the parameter range explored, is found to be 
always subrelativistic, owing to the large neutrino-driven mass flux expelled along the 
field lines. The slow magnetosonic point in this case is located very near the surface, 



T =2 MeV T =4 MeV 




z/R z/R fl 

FIGURE 1. Profiles of various quantities in the sub-slow magnetosonic region, computed using a split 
monopole magnetic field with tan 9 = 0.2, where 9 is the inclination angle of the field line with respect to 
the symmetry axis, Keplerian rotation Q. = Q.^, and surface magnetic field B p q = — lOB^o = 10 15 G. Each 
panel corresponds to a run with a different neutrino luminosity (indicated in terms of the effective black- 
body temperature T v ). The quantities plotted in each panel are: the slow magnetosonic Mach number 
M$m (solid line), the dimensionless entropy per baryon s (dashed line), the temperature T in units of the 
initial temperature Tq (dotted-dashed line), and the energy deposition per baryon per unit length along the 
streamline measured in units of m p c 2 /Rq, dS jdl (dotted line). All quantities are given as functions of the 
normalized height above the disk midplane z/Rq- 
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FIGURE 2. Dependence of mass-to-magnetic flux ratio Tj (left axis) on the inclination angle 9 of the 
split-monopole field line. The corresponding values of the mass flux M are indicated on the right axis 



at z/Rq « 1. A systematic study of the dependence of r\ on the various parameters is 
presented in Ref. [9]. 



CONCLUSION 



Our main conclusion is that ejection of relativistic outflows from the innermost disk 
radii, within several r s or so, is possible in principle for certain magnetic field con- 
figurations even in non-rotating black holes, provided the neutrino luminosity emitted 
by the disk does not exceed 10 52 ergs s _1 or so, and the magnetic field is sufficiently 
strong (B p ~ 10 15 G) . The neutrino-driven mass flux depends rather sensitively on the 
neutrino luminosity, on magnetic field geometry in the sub-slow magnetosonic region, 
and on the angular velocity of magnetic flux surfaces, but is highly insensitive to the 
strength of poloidal and toroidal magnetic field components near the surface provided 
the Alfven Mach number is sufficiently small there. The picture often envisioned, of an 
ultra-relativistic core surrounded by a slower, baryon rich wind J2,[l5|] seems a natural 
consequence of a neutrino-assisted MHD disk outflow. 

The hyperaccretion process is likely to be intermittent, leading to temporal changes 
in the neutrino luminosity. This should result in large variations of the Lorentz factors of 
consecutive fluid shells expelled from the disk in the polar region, owing to the sensitive 
dependence of mass loading on L v . In this situation we expect strong shocks to form in 
the outflow. If the polar disk outflow is associated with the GRB-producing jet, then the 
observed gamma-ray emission can be quite efficiently produced behind those shocks. 

Finally, the conditions we find to be optimal for launching an ultra-relativistic jet in 
the polar region, are also the conditions favorable for large neutron-to-proton ratio in 
the disk. However, for the steady flow considered above we estimate the ratio of the 
neutronization timescale and flow time to be of order unity, and so the electron fraction 
is expected to evolve as the flow accelerates. Detailed analysis of the composition profile 
in the outflow is left for future work. 
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